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Abstract: Recently, refinements have been made on both the theoretical and experimental 
determinations of i.) the mass of the lightest Higgs scalar (m/J, ii.) the relic density of 
cold dark matter in the universe (QcDMh 2 ), Hi.) the branching fraction for radiative B 
decay BF(b — > sj), iv.) the muon anomalous magnetic moment {an), and v.) the flavor 
violating decay B s — > Each of these quantities can be predicted in the MSSM, 

and each depends in a non-trivial way on the spectra of SUSY particles. In this paper, 
we present updated constraints from each of these quantities on the minimal supergravity 
(mSUGRA) model as embedded in the computer program ISAJET. The combination of 
constraints points to certain favored regions of model parameter space where collider and 
non-accelerator SUSY searches may be more focussed. 
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1. Introduction 

The search for weak scale supersymmetric matter is one of the prime objectives of present 
and future collider experiments. Particle physics models including supersymmetry solve a 
host of problems occurring in non-supersymmetric theories, and predict a variety of new 
matter states — the sparticles — at or around the TeV scale 0. Supersymmetric models can 
be classified by the mechanism for communicating SUSY breaking from the hidden sector to 
the observable sector. Possibilities include gravity mediated SUSY breaking (SUGRA)|2|, 
gauge mediated SUSY breaking (GMSB)[|, anomaly mediated SUSY breaking (AMSB)f§ 
and gaugino mediated SUSY breaking (inoMSB)|5[. Of these, the SUGRA models may 
be perceived as the most conservative, since they do not require the introduction of either 
extra dimensions or new messenger fields, and because gravity exists. 

The so-called minimal supergravity (mSUGRA) model (sometimes also referred to 
as the CMSSM model) has traditionally been the most popular choice for phenomeno- 
logical SUSY analyses. In mSUGRA, it is assumed that the Minimal Supersymmetric 
Standard Model (MSSM) is valid from the weak scale all the way up to the GUT scale 
Mqut — 2 x 10 16 GeV, where the gauge couplings g\ and gi unify. In many of the early 
SUGRA models^, a simple choice of Kahler metric G\ and gauge kinetic function f ab led 
to universal soft SUSY breaking scalar masses (mo), gaugino masses (m^j) and A-terms 
(^4o) a t Mgut- This assumption of universality in the scalar sector leads to the phe- 
nomenologically required suppression of flavor violating processes that are supersymmetric 
in origin. However, there is no known physical principle which gives rise to the desired 
form of G\ and Jab] indeed, for general forms of G\ and Jab, non- universal masses are 
expected ||]. In addition, even if nature did select a SUGRA model leading to tree level 
universality, quantum corrections would (without further assumptions) lead to large devi- 
ations from universality Q. Hence, the universality assumption nowadays is regarded as 
being ad hoc — entirely motivated by the phenomenological need for suppression of flavor 
violating processes in the MSSM. 

In the mSUGRA model, we thus assume universal scalar masses, gaugino masses and 
A-terms. We will also require that electroweak symmetry is broken radiatively (REWSB), 
allowing us to fix the magnitude, but not the sign, of the superpotential Higgs mass term \i 
so as to obtain the correct value of Mz- Finally, we trade the bilinear soft supersymmetry 
breaking (SSB) parameter B for tan f3 (the ratio of Higgs field vacuum expectation values). 
Thus, the parameter set 

mo, rni/2, Aq, tan/3, and sign(fi) (1-1) 

completely determines the spectrum of supersymmetric matter and Higgs fields. 

In our calculations, we use the program ISASUGRA to calculate the SUSY particle 
mass spectrum. ISASUGRA is part of the IS A JET || package. Working in the DR regu- 
larization scheme, the weak scale values of gauge and third generation Yukawa couplings 
are evolved via 2-loop RGEs to Mgut- At Mgut, universal SSB boundary conditions 
are imposed, and all SSB masses along with gauge and Yukawa couplings are evolved to 
the weak scale M wea k- Using an optimized scale choice Qsusy = ^ m t L m t R ' the RG- 
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improved one-loop effective potential is minimized and the entire spectrum of SUSY and 
Higgs particles is calculated. Our values of m^ are in close accord with those generated 



by the FeynHiggsFast program^. Yukawa couplings are updated[10| to account for SUSY 
threshold corrections, and the entire parameter set is iteratively run between M„ ea i. and 
Mqut until a stable solution (within tolerances) is obtained. 

Once the SUSY and Higgs masses and mixings are known, then a host of observables 
may be calculated, and compared against experimental measurements. The most important 
of these include: 

• lower limits on sparticle and Higgs boson masses from new particle searches at LEP2, 

• the relic density of neutralinos originating from the Big Bang, 

• the branching fraction of the flavor changing decay b —* 57, 

• the value of muon anomalous magnetic moment = ^ 2 2 ~ >M and 

• the lower bound on the rate for the rare decay B s — * 

Our goal is to delineate the mSUGRA parameter space region consistent with all these 
constraints. Similar studies have recently been presented in Refs. 12, 13, 14, |i"5| . In 



this report we examine these constraints using the updated ISASUGRA package, which 
is convenient for explicit event generation for various colliders using the ISAJET v7.63 
program, which includes several improvements over previous versions. The most important 
of these improvements are the evaluation of the bottom Yukawa coupling and the value of 
rriA, especially for large values of tan/5. Complete 1-loop self energy corrections as given 
in Ref. fl(| have been incorporated. 1 Comparisons of ISAJET v7.63 with other similar 
codes are available in Ref. [jig] . 

Returning to our analysis, we also incorporate a new calculation of the neutralino relic 
density ^^Ji? that has recently become available [17]. In Ref. [17], all relevant neutralino 



annihilation and co- annihilation processes are calculated, and the neutralino relic density 



is evaluated using relativistic thermal averaging (see also Refs. [ 18 ] and recently Belanger 
et al, Ref. gjj). The latt er is especially important in evaluating the relic density when 
s-channel annihilation resonances occur, as in the mSUGRA model at large tan (3, when 
Z1Z1 — > A, H — > //, where the fs are SM fermions[p^, [2C| , jHj , [22^ , 23] . We also present 
improved b — > S7 branching fraction predictions in accord with the current ISAJET release. 
We also discuss constraints imposed by the measurement of the muon anomalous magnetic 
moment, using updated calculations from two somewhat different analyses. Finally, we 
delineate the region of mSUGRA parameter space excluded by the CDF lower limit |p4|| on 
the branching fraction of B s — * ■ This constraint is important only for very large 

values of tan/3 pq] . 

Within the mSUGRA framework, the parameters mo and m^ 2 are the most impor- 
tant for fixing the scale of sparticle masses. The mo — mi/2 plane (for fixed values of other 



1 On the technical side, the numerical precision in ISAJET has been improved to facilitate a better 
analysis near the boundary of the region excluded by electroweak symmetry breaking constraints. 
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parameters) is convenient for a simultaneous display of these constraints, and hence, of pa- 
rameter regions in accord with all experimental data. Physicists interested in the mSUGRA 
model may wish to focus their attention on these regions. We also present five mSUGRA 
model cases illustrating distinctive characteristics of the mSUGRA particle spectrum for 
parameter choices which are consistent with all experimental constraints. 

The remainder of this paper is organized as follows. In Sec. 0, we discuss the various 
constraints on the mSUGRA model, and present some details of our calculations. In Sec. |3|, 
we show our main results as regions of the tjiq vs. parameter space plane for different 
values of tan/3 and sign of fi. Our conclusions and sample points are presented in Sec. |I[ 



2. Constraints and calculations in the mSUGRA model 
2.1 Constraints from LEP2 searches 

The LEP2 collaborations have finished taking data, and significant numbers of events were 
recorded at e + e~ CM energies ranging up to y^s ~ 208 GeV. Based on negative searches 
for superpartners at LEP2, we require 



m 



w > 103.5 GeVfH and 



• mg L R > 99 GeV provided — < 10 GeV[27], which is the most stringent of 
the slepton mass limits. 

The LEP2 experiments also searched for the SM Higgs boson. In addition to find- 
ing several compelling signal candidates consistent with ~ 115 GeV, they set a limit 
m H SM > H4.1 GeV|28|]. In our mSUGRA parameter space scans, the lightest SUSY Higgs 
boson h is almost always SM-like. The exception occurs when the value of ttia becomes 
low, less than 100 — 150 GeV. Then there exists a near mass degeneracy between h and H, 
and the SM Higgs is a mixture of these. This case arises at very large values of tan (3. For 
clarity, we show contours where 

• m h > 114.1 GeV, 

and will direct the reader's attention to any regions where this bound might fail. 
2.2 Neutralino relic density 

Measurements of galactic rotation curves, binding of galactic clusters, and the large scale 
structure of the universe all point to the need for significant amounts of cold dark matter 
(CDM) in the universe. In addition, recent measurements of the power structure of the 
cosmic microwave background, and measurements of distant supernovae, point to a cold 



dark matter density [29] 



• 0.1 < tt C DMh z < 0.3 . 

The lightest neutralino of mSUGRA is an excellent candidate for relic CDM particles in 
the universe. The upper limit above represents a true constraint, while the corresponding 
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lower limit is flexible, since there may be additional sources of CDM such as axions, or 
states associated with the hidden sector and/or extra dimensions. 

To estimate the relic density of neutralinos in the mSUGRA model, we use the recent 



calculation in Ref. jlT]] . In Ref . 17 ] , all relevant neutralino annihilation and co-annihilation 
reactions are evaluated at tree level using the CompHEP|}0) program for automatic eval- 
uation of the associated 7618 Feynman diagrams. The annihilation cross section times 



velocity is relativistically thermally averaged[31], which is important for obtaining the cor 



rect neutralino relic density in the vicinity of annihilations through s-channel resonances. 
2.3 The b — > sj branching fraction 



The branching fraction BF(b — > S7) has recently been measured by the BELLE [32], 
CLEO|33| and ALEPHj34| collaborations. Combining statistical and systematic errors 
in quadrature, these measurements give (3.36 ±0.67) x 1(T 4 (BELLE), (3.21 ±0.51) x 10" 4 
(CLEO) and (3.11 ± 1.07) x 10~ 4 (ALEPH). A weighted averaging of these results yields 
BF(b -» 87) = (3.25 ± 0.37) x 10" 4 . The 95% CL range corresponds to ±2cx away from 
the mean. To this we should add uncertainty in the theoretical evaluation, which within 
the SM dominantly comes from the scale uncertainty, and is about 10%. 2 Together, these 
imply the bounds, 



2.16 x 10" 4 < BF(b -► S7) < 4.34 x 10" 4 



Other computations of the range of BF{b — > sj) include for instance Ellis et al. ]12|: 
2.33 x 10- 4 < BF{b -> 57) < 4.15 x 10~ 4 , and Djouadi et a/.@: 2.0 x 10~ 4 < BF(b -» 
57) < 5.0 x 10 -4 . In our study, we simply show contours of BF(b — ► S7) of 2, 3, 4 and 
5 x 10 -4 , allowing the reader the flexibility of their own interpretation. 

The calculation of BF(b — ► 57) used here is based upon the program of Ref. p6[ . 
That calculation uses an effective field theory approach to evaluating radiative corrections 
to the b — > sj decay rate. In running from Mqut to M wea k, when any sparticle threshold 
is crossed, the corresponding sparticle is integrated out of the theory, and a new basis 
of decay-mediating operators multiplied by Wilson coefficients (WCs) is induced. The 



evolution of the WCs can be calculated by RG methods. We adopt the Anlauf procedure^] 
in our calculation, which implements a tower of effective field theories, corresponding to 
each sparticle threshold which is crossed. This procedure sums large logarithms that can 
occur from a disparity between different scales involved in the loop calculations. In our 
calculations, we implement the running 6-quark mass including SUSY threshold corrections 
as calculated in ISASUGRA; these effects can be important at large values of the parameter 



tan/^J^, 39]. Once the relevant operators and Wilson coefficients are known at Q = Mw, 
then the SM WCs are evolved down to Q = nib via NLO RG running. At nib, the BF(b — ► 
57) is evaluated at NLO, including bremsstrahlung effects. Our value of the SM b — > S7 
branching fraction yields 3.4 x 10 -4 , with a scale uncertainty of 10% 



0. 



2 We caution the reader that the SUSY contribution may have a larger theoretical uncertainty, particu- 
larly if tan p is large. An additional theoretical uncertainty that may increase the branching ratio in the 
SM is pointed out in Ref. EH). 
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2.4 Muon anomalous magnetic moment 

The muon anomalous magnetic moment = ^ ^ has been recently measured to high 
precision by the E821 experiment Q: a M = 11659202(14)(6) x 1(T 10 . In addition, addi- 
tional data analyses should soon be finished, and we may anticipate a further reduction in 
the experimental error by a factor of 2. The initially reported^] 2.6<r deviation from the 



SM value of Ref.[41] has since been tempered somewhat by correcting the sign of the SM 
light-by-light contribution to q» [^| . A correction of the sign of the hadronic light by light 
contribution reduces the significance of the deviation from the SM to 1.6a, i.e. at 2a: 



-6 < 6a^ x 10 10 < 58 (CM). 



An alternative evaluation of theory uncertainties in the SM calculation by Melnikov[43] 
leads to (including the LBL correction): 

• -29.9 < 5a^ x 10 10 < 62.3 (Melnikov). 3 

In view of the theoretical uncertainty, and the impending new experimental analysis, we 
only present contours of Sa^, as calculated using the program developed in p5| , and leave 
it to the reader to decide the extent of the parameter region allowed by the data. 

2.5 B s — > fi + [i~ decay 

While all SUSY models contain two doublets of Higgs superfields, there are no tree level 
flavor changing neutral currents because one doublet H u couples only to T3 = 1/2 fermions, 
while the other doublet couples just to T3 = —1/2 fermions. At one loop, however, 
couplings of H u to down type fermions are induced. These induced couplings grow with 
tan p. As a result, down quark Yukawa interactions and down type quark mass matrices 
are no longer diagonalized by the same transformation, and flavor violating couplings of 
neutral Higgs scalars h, H and A emerge. Of course, in the limit of large dia, the Higgs 
sector becomes equivalent to the SM Higgs sector with the light Higgs boson h = Hsm, and 
the flavor violation decouples. The interesting thing is that while this decoupling occurs 
as rriA — > 00, there is no decoupling for sparticle masses becoming large. 

An important consequence of this coupling is the possibility of the decay B s — > 
whose branching fraction has been experimentally bounded by CDF[[24|] to be: 

• BF{B S fi+fj,-) < 2.6 x 1(T 6 , 

mediated by the neutral states in the Higgs sector of supersymmetric models. While this 
branching fraction is very small within the SM (BF SM {B S -> fi + fi~) - 3.4 x 1(T 9 ), the 
amplitude for the Higgs-mediated decay of B s grows as tan 3 (3 within the SUSY framework, 



and hence can completely dominate the SM contribution if tan (5 is large. Several groups [25] 



have analyzed the implications of this decay within the mSUGRA framework. A subset of 



3 Melnikov — who uses the analysis of Ref.|44[ which does not use tau decay data for the evaluation of 
the hadronic vacuum polarization and has a more conservative error on the light by light contribution — 
finds <5a M x 10 10 = 16.2 ± 14.0| s t a t ± 6.0| 3ys ± 15.6|t/ieorj/- The conservative "2<r" range reported here can 
be obtained by linearly combining the theory error with the 2<r experimental error. 
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us [46 1 have recently performed an independent analysis of this decay. In the following, we 
use the results of this analysis to delineate the region of mSUGRA parameters excluded 
by the CDF upper limit on its branching fraction. 

Tevatron experiments should be able to probe this decay in the near future. With an 
integrated sample of 2 fb _1 they should be sensitive to a branching fraction for B s — ► 
down to ~ 1CP 7 . With a still bigger data sample (that is expected to accumulate before 
the Large Hadron Collider begins operation) the sensitivity should be even greater. 

3. Results 

Our first results are plotted in Fig. [l|. Here, we show the ttiq vs. mi/2 plane for Aq = 0, 
tan (3 = 10 and /u < 0. The red shaded regions are excluded either due to a lack of REWSB 
(right-hand side), or a stau LSP (left-hand side). The magenta region is excluded by 
searches for charginos and sleptons at LEP2. The region below the red contour is excluded 
by LEP2 Higgs searches, since here rrih < 114.1 GeV. In addition, we show regions of 
neutralino relic density with 0.1 < £lg h < 0.3 (green), which is favored by cosmological 
observations, and also 0.3 <VL^Ji 2 <1 (blue) and 0.02 < Qg^h 2 < 0.1 (yellow). The bulk 
of the unshaded region in the center of the plot has Q,^ h 2 > 1, and would thus be excluded 
since the age of the universe would be less than 10 billion years. The magenta contours 
denote values of BF(b — » sj) = 4 and 5 x 10 -4 . Finally, the blue contours denote values 
of 5a^ = —30, —10, —5 and —2 x 10~ 10 , moving from lower left to upper right. There is no 
constraint arising from B s —* decay at tan (3 = 10. 

An intriguing feature of the plot is that the lower left green shaded region of desirable 
relic density, where neutralinos mainly annihilate via t-channel slepton exchange to lepton- 
anti-lepton pairs is essentially excluded by the m^, b — > s~f and 5a^ constraints. That leaves 
two allowed regions with a preferred relic density: one runs near the stau LSP region, where 
fi — Z\ co-annihilation effects reduce an otherwise large relic density (as pointed out by 
Ellis et aZ.[^7|). This region has a highly fine-tuned relic density, since a slight change in 



rriQ leads to either too light or too heavy of a fi mass to give 0.1 < Sl/i 2 < 0.3[48, 17]. 
The other region runs parallel to the REWSB excluded region, and occurs when the Z\ 
has a sufficiently large higgsino component that annihilation into WW, ZZ and Zh pairs 



reduces the relic density 49, 17]. This region corresponds to what is known as "focus point" 
SUSY, and since tuq is large, SUSY scalar masses are also large, leading to some degree 
of suppression of FCNC and CP violating processes [50]. The narrowness of the region 



indicates again that some fine-tuning of parameters is needed to achieve the right relic 
density, although the amount of fine-tuning is less than in the f\ co-annihilation case. 

A similar plot is shown in Fig. ^, but in this case for [i > 0. Much of the labeling is 
similar to Fig. Q, although now the b — > sj contour denotes a branching fraction of 3 x 10 -4 . 
In this case almost the entire plane shown is in accord with the measured branching fraction 
for this decay. In addition, the blue contours denote values of 5a^ = 60, 40, 20, 10, 5 and 
2 x 10~ 10 . Constraints from 5a^ as well as from B s — > fi + \i~ are not relevant for this case. 

For tan/3 = 10 and /i > 0, the slepton annihilation region of relic density has a small 
surviving region just beyond the Higgs mass contour. For the most part, to attain a 
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tanfi = 10, A = 0, [L<0 




1000 1500 2000 2500 3000 
m Q ( GeV) 



Figure 1: Plot of constraints for the mSUGRA model in the mo vs. m 1 / 2 plane for tan/? = 10, 
Aq = and /i < 0. We show regions of CDM relic density, plus contours of m/j = 114.1 GeV, 
contours of muon anomalous magnetic moment (xlO 10 ) and contours of b — > S7 branching 
fraction (xlO 4 ). 



preferred value of neutralino relic density, one must again live in the stau co-annihilation 
region, or the focus point region. A final possibility is to be in the slepton annihilation 
region, but then the value of should be slightly beyond the LEP2 limit; in this case, a 
Higgs boson signal may be detected in Run 2 of the Fermilab Tevatron[^]. 

We next turn to the rriQ vs. m 1 / 2 plane for tan/3 = 30 and /i < 0. The gray region 
in the bottom left corner of the plot is excluded because rn~^ < 0. In this case, the area 
of the green region of relic density in the lower-left has expanded considerably owing to 
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tanfi = 10, A = 0, jo, > 




m n ( GeV) 



Figure 2: Same as Fig. [j], but for fj, > 0. 

enhanced neutralino annihilation to bb and rf at large tan (3. Both lighter values of m^ x 
and and also large r and b Yukawa couplings at large tan (3 enhance these t-channel 
annihilation rates through virtual staus and sbottoms. Unfortunately, the region excluded 
by BF(b — > sj) and by Sa^ (even with the conservative constraint of Ref.f4|[]) also expands, 
and most of the cosmologically preferred region is again ruled out. As before, we are left 
with the corridors of stau co-annihilation and the focus point scenario as the only surviving 
regions. 

The corresponding plot is shown for tan/3 = 30 but fi > in Figure |j. In this case, 
the magenta contours of BF{b — > 57) correspond to 2 and 3 x 10 -4 . Thus, the lower left 
region is excluded since it leads to too low a value of BF(b — > 57). The 8a^ contours begin 
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tanfi = 30, A = 0, [L<0 




m n ( GeV) 



Figure 3: Same as Fig. [j], but for tan/3 = 30 and /i < 0. The light blue contour labeled 0.1 denotes 
where B(B S — ► = 0.1 x 10~ 7 . In subsequent figures these branching fractions contours are 

all labeled in units of 10 -7 . 

from lower left with 60 x 10~ 10 , then proceed to 40, 20, 10, 5 and 2 x 10~ 10 . A fraction of 
the slepton annihilation region of relic density is excluded also by too large a value of 5a M . 
Of course, a reasonable relic density may also be achieved in the stau co-annihilation and 
focus point regions of parameter space. 

Next, we examine the mSUGRA parameter plane for very large values of tan = 45 and 
[i < 0. If we take tan (5 much bigger than 45 for this sign of fi, the entire parameter space is 
excluded due to lack of REWSB. The gray and red regions are as in previous figures. The 
blue region is excluded because m 2 A < 0, denoting again a lack of appropriate REWSB. 
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tanfi = 30, A = 0, jo, > 




m n ( GeV) 



Figure 4: Same as Fig. |], but for tan f3 = 30 and [i > 0. 

The inner and outer red dashed lines are contours of ttia = 100 and ttia = 200 GeV, 
respectively. The former is roughly the lower bound on m\ from LEP experiments. In 
between these contours, h is not quite SM-like, and the mass bound from LEP may be 
somewhat lower than = 114.1 GeV shown by the solid red contour, but outside the 
200 GeV contour this bound should be valid. 

The tiny black region in the figure is excluded by constraints |52jj on the width of the 
Z boson (specifically, the decay Z — > hA + HA leads to too large a value for Tz)- We also 
see that much of the lower-left region is excluded by too high a value of BF(b — * sj) and 
too low a value of 5a In addition, in this plane, the experimental limit on B s — * 
enters the lower-left, where values exceeding 26 x 10~ 7 are obtained. It seems that in the 
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upper region which is favored by the b — > sj constraint, detection of B s — > at the 

Tevatron will be quite challenging. 

In this figure, the relic density regions are qualitatively different from the lower tan (3 
plots. A long diagonal strip running from lower-left to upper-right occurs because in this 
region, neutralinos annihilate very efficiently through s-channel A and H Higgs graphs, 
where the total Higgs widths are very large due to the large b and r Yukawa couplings for 
the high value of tan/3 in this plot. Adjacent to this region are yellow and green shaded 
regions where neutralino annihilation is still dominated by the s-channel Higgs graphs, but 
in this case the annihilation is somewhat off-resonance. The A and H widths are so large, 
typically 20 — 60 GeV, that even if \2m^ — m^ H ^\ is relatively large, efficient annihilation 
can still take place. The relic density changes so slowly on the flanks of the annihilation 
corridor that little fine tuning of parameters is needed to achieve a favored value of tig h 2 . 

For the case of fi > 0, tan (3 values ranging up to 60 can be allowed, although the 
mSUGRA parameter space becomes very limited for tan/3 ~ 55. Hence, we show in Fig. 
H the mSUGRA parameter space plane for tan /? = 52 and /i > 0. In this parameter 
plane, the relic density annihilation corridor occurs near the boundary of the excluded 
f\ LSP region. The width of the A and H Higgs scalars is very wide, ranging from 30 
GeV for m 1/2 ~ 400 GeV, to 110-130 GeV for m 1/2 ~ 2000 GeV. Efficient s-channel 
annihilation through the Higgs poles can occur throughout much of the upper region of 
allowed parameter space. But the annihilation is not overly efficient due to the extreme 
breadth of the Higgs resonances. In fact, none of the entire parameter plane is excluded 
by h 2 > 1 ! The change in h 2 is so slow over almost the entire parameter plane that 
little fine-tuning occurs. 

In much of the region with m^ 2 < 400 GeV, the value of BF(b — > sj) is below 
2 x 10 -4 , so that some of the lower green relic density region where annihilation occurs 
through i-channel stau exchange is excluded. In contrast, the value of 5a^ is in the range 
of 10 - 40 x 10~ 10 , which is in accord with the E821 measurement. The value of mh is 
almost always above 114.1 GeV, and the BF(B S — ► n + is always below 10~ 7 , and could 
(if at all) be detected with several years of main injector operation. 

Aside from the somewhat low value of BF(b — * sj) in the lower left, much of this 
plane represents a very attractive area of mSUGRA model parameter space. If the model 
parameters are indeed in this range, the Tevatron signal for B s — > ji + [i~ will be small, and 
BF(b — > S7) will turn out somewhat below the SM prediction, while 5a^ will be somewhat 
above the SM value. 

4. Conclusions 

In this paper, we have presented updated constraints on the mSUGRA model from i.) the 
LEP2 constraints on sparticle and Higgs boson masses, ii.) the neutralino relic density 
Qg^h 2 , Hi.) the branching fraction BF(b — > sj), iv.) the muon anomalous magnetic 
moment and v.) the leptonic decay B s — > /i + /x~. Putting all five constraints together, 
we find favored regions of parameter space which may be categorized by the mechanism 
for annihilating relic neutralinos in the early universe: 
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tanfi = 45, A Q = 0, [L<0 




500 1000 1500 2000 2500 3000 3500 4000 



m n ( GeV) 



Figure 5: Same as Fig. but for tan j3 = 45 and \i < 0. The inner and outer red dashed lines are 
contours of tua = 100 and rriA = 200 GeV, respectively. 

• 1. annihilation through t-channel slepton exchange (low mo and m^j), 

• 2. the stau co-annihilation region (very low mo but large m^ 2 ), 

• 3. the focus point region (large mo but low to intermediate m^) and 

• 4. the flanks of the neutralino s-channel annihilation via A and H corridor at large 
tan (3 when Ta and Yh are very large. 

In previous years, there may have been a preference for region 1. as offering the most 
natural channel for obtaining a reasonable value of relic density. However, recently much of 
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tanfi = 52, A Q = 0, [L> 




a. susy 



- m h =114.1GeV 
0.3 <Q.h 2 <1.0 
0.1 <Q.h 2 <0.3 
0.02 <Q.h 2 <0.1 

Br(b^>sy) 



1000 1500 2000 2500 3000 3500 4000 
m Q ( GeV) 



Figure 6: Same as Fig. |], but for tan (3 = 52 and [i > 0. 



this region is ruled out by a combination of LEP2 limits on at low tan /3, too high (for 
fi < 0) or too low (for \i > and large tan 0) a value of BF(b — > S7), and too low a value 
of a M (for fj, < and intermediate to large tan/3). In addition, the CDF BF(B S — > fi + 
constraint is starting to become important for /x < and large tan/3. Nevertheless, some 
of region 1 remains viable, especially for \i > 0, where we expect a lightest Higgs boson 
just beyond the bounds from LEP2. 



The stau co-annihilation region[47] 2. is intriguing because one can always take m^j 
large enough for any tan (3 value to evade constraints on deviations from SM predictions. 
However, this region is exceptionally narrow in the parameter mo, and slight deviations 
cause either too high or too low a value of relic density. This indicates a high degree of 
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fine-tuning in the determination of the relic density in this region. 



The focus point region[5C] 3. also leads to a reasonable relic density, this time because 
the higgsino component of Z\ is large enough that efficient annihilation can occur to WW, 
ZZ and ZH states. It may also be preferable based on possibly low electroweak fine-tuning, 
and because matter scalar masses are high enough to offer some degree of suppression of 
SUSY induced FC and CP violating processes 5C]. This region also suffers some degree of 
fine-tuning of the relic density, since too high or too low a higgsino component of Z\ can 
result in too low or too high a value of relic density. For tan (3 ~ 10 and > 0, the focus 
point region is in accord with all constraints. For \x < 0, the focus point region usually 
gives too high a value of BF{b — > s-y); for ^ > 0, we would expect ultimately experiment to 
measure a somewhat lower value of BF(b — > sy) than the SM prediction, and a somewhat 
higher value of . 

Finally, at very large tan (3 ~ 45 — 55 there can exist wide regions of parameter space 
where Z\Z\ annihilation can occur in the early universe through very broad s-channel 
A and H resonances, giving rise to a reasonable value of relic density |20|, [Eg]: region 4.. 
Unfortunately, much of this region is excluded for fi < by too large a value of BF(b — > sy). 
Here also a rather low a value of less than the SM prediction is generated. In this case, 
very large values of mo and m^/2 are needed to escape experimental constraints, perhaps 
placing the SUSY spectrum in conflict with naturalness bounds [53]. 

For very large tan (3 and /j > 0, however, broad regions of parameter space can be 
found with a reasonable relic density, and also which are in accord with all other low energy 
constraints. In this region, we expect BF(b — > 57) somewhat below the SM prediction, and 
somewhat above the SM prediction. The current BF(b — ► 57) measurement suggest 
m i/2 values ~ 500 GeV, giving rise to sparticles typically at the TeV scale or beyond. 

To summarize these various regions, we present in Table |l] five parameter space points 
indicating the SUSY spectrum that might occur in each region. We also list the relic 
density, BF(b — > 57), a M and BF(B S — » values. 

Point 1 occurs in region 1., and is characterized by light sparticle masses, especially tau 
sleptons. The light Higgs scalar h is slightly beyond the LEP2 bound. The search channel 
at the Fermilab Tevatron would be pp — ► W1Z2X, with Z2 — ► tt\ and W\ — ► fu T though 
detection appears to be difficult [54]. The CERN LHC would be awash in signals [55], and 
many new SUSY states would be accessible to a linear collider (LC) with ^/s ~ 500 — 800 
GeV|6|. 

Point 2 is in the stau co-annihilation region 2., but at low enough mi/2 that i-channel 
stau exchange is still important in the Z%Z\ annihilations. It is characterized by a rather 
small mass gap between fi and Z\. Only the Higgs boson h would be observable at the 
Tevatron. At the LHC, a variety of leptonic signatures would occur in gluino and squark 
cascade decays, including a high rate of r production. The LC would have to operate 
slightly above yfs ~ 500 GeV to access even the first SUSY states. 

Point 3 in the focus point region 3. has TeV scale scalars, but light W\ and Zi- The 
LSP is a mixture of gaugino-higgsino. The W\ decays into three bodies dominated by W* 
exchange, and the Z2 decay is dominated by Z* exchange. Because this point is in the 
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region with small \fj,\ (but with a not too small — ti.^), a low rate of trilepton events 
may be accessible to Fermilab Tevatron experiments. At the CERN LHC, gg and W\Zi 
production would be dominant. A LC would be able to find most of the charginos and 
neutralinos, since these are the lowest lying states. 

Point 4a lies in region 4. on the edge of the Higgs annihilation corridor for large tan j3 
and \x < 0. To evade constraints, the SUSY spectrum is very heavy: only the Higgs h 
would be seen at the Tevatron or a LC, and in fact it would even be challenging to discover 
SUSY at the LHC. This point is disfavored by naturalness arguments. 

Point 4b also in region 4. lies in the Higgs annihilation corridor for large tan (3 and 
[i > 0. A significantly lighter spectrum can be tolerated than in the 4a case, although 
SUSY scalars are still in the TeV range. Only the h would be accessible at the Tevatron. 
SUSY signals corresponding to point 4b should be readily visible at the LHC, although a 
LC with y/s > 1 TeV would be required to access even the lowest lying SUSY states. 

To summarize, we find the five constraints considered in this paper to be highly re- 
strictive. Together, they rule out large regions of parameter space of the mSUGRA model, 
including much of the region where i-channel slepton annihilation of neutralinos occurs in 
the early universe. The surviving regions 1.-4. have distinct characteristics of their SUSY 
spectrum, and should lead to distinct SUSY signatures at colliders. 
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Table 1: Representative weak scale sparticle masses (in GeV units) and parameters for five selected 
mSUGRA models. We use Aq = and m t = 175 GeV. The value of \x is also shown since it is 
sometimes regarded as a measure of fine tuning. 



parameter 






value 






•point 
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4a 


4b 




100 


165 


1200 


2750 


800 


^1/2 


300 


550 


250 


1800 


800 


tan(/3) 


10 


10 


10 


45 


52 


sgn(n) 


1 


1 


1 


-1 
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m- g 


701.4 


1225.6 


658.0 


3810.8 


1757.6 


mu L 


630.7 


1099.1 


1271.1 


4185.0 


1715.3 


m u R 


611.1 


1060.0 


1269.1 


4080.7 


1662.4 


TOr 


635.6 
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m ~ bl 
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1484.4 
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1539.4 


m h 
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858.5 


825.2 
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1328.2 


m h 
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3529.6 
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on c /i 

395.4 


i i no n 

1198.0 


2750.8 


834.0 


m fl 


1 A T EC 

147.5 


257.6 


1191. (J 
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260.9 
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BF(b -► 87) x 10 4 


3.12 


3.46 


3.20 


3.92 


2.85 


a s ^ USY x 10 10 
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